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f  preparation  of  this  report. 
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A  TRANS- IONOSPHERIC-SIGNAL  SPECIFICATION 
FOR  SATELLITE  C3  APPLICATIONS 


I .  INTRODUCTION 

Proper  design  of  radio  frequency  systems  that  must  operate 
through  structured  or  varying  plasmas  require  an  accurate  and  easily 
iraplementabte  signal  specification.  The  specification  should  be  a 
"reasonable  worst  case"  in  that  any  system  which  operates  with  the 
specified  signal  can  also  operate  under  any  likely  disturbed  condition. 

The  following  specification  is  the  result  of  a  dedicated 

program  sponsored  by  the  Defense  Nuclear  Agency  and  the  Air  Tore© 

Weapons  Laboratory.  Other  taaj^r  contributors  were  the  Naval  Research 

Laboratory  and  the  Air  Force  Geophysics  Laboratory.  The  program  included 

investigations  in  the  mechanisms  that  eause  structured  plasmas,  in  the 

propagation  of  electromagnetic  signals  through  structured  plasmas,  and 

in  the  effects  of  scintillated  or  otherwise  distorted  signals  on  typical 
5 

satellite  C  links.  Field  experiments  and  measurement  programs  were 
executed  to  verify  theoretical  plasma  structure  and  radio  propagation 
predict  ions  and  to  characterise  the  morphology  of  the  natural  ionosphere. 
Implicit  in  the  resulting  specification  are  several  considerations  and 
assumptions.  Rather  than  discuss  all  of  the  supporting  logic  leading  to 
this  specification,  a  selected  bibliography  is  included.  ether  supporting 
iaforaatioa  is  available  in  the  classified  literature.  It  is  ftsted, 
however,  that  this  specification  should  be  adequate  for  all  satellite  C* 
and  data  link  applications. 


? 

I 


} 
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2. 


SIGNAL  EFFECTS  FORMALISM 


The  effects  of  a  disturbed  ionospheric  channel  are  represented 
by  the  channel  impulse  response  function. 


U(t) 


oo 

J dt  h(t,t)  S(t-t) 


(1) 


S(t)  =  transmitted  signal 

R(t)  =  received  signal  (complex  number  representation) 
h(t,r)  *  chaimel  impulse  response  function 


Tl*e  specification  channel  impulse  response  function  is 


shore  f  »  carrier  frequency  (hs) 

ft 

e  a  light  speed  (Sxl©  a/sec) 

sit)  »  propagation  path  length  (a) 

?  9  classical  electron  radius  Cl.S-xlQ  a) 

a  % 

Mft)  9  total  electron  content  ia""l 


far  ather  than  circular  polarisation,  fustian  1  should  be  applied  to  each 
polarisation  state  and  N(t)  used  to  calculate  the  Faraday  rotation  effects. 
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(3) 


A(t)  *  exp( 


•o.nsKA)i 


exp£-e2/2o2) 


*  absorption  (db) 

(1(9)  *  antenna  gain  function 

2  2 
ea  *  energy  angle  of  arrival  variance  (rad  ) 

tj 


S» 


(va  tm 
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|df  h  £f,T)exp(i2wAfT*i2*ft) 


(4) 


Equation  3  assumes  that  the  transmitter  is  on  the  antenna  axis  and  that 
the  antenna  pattern  is  cylindrical ly  symmetric  about  that  axis.  If  the 
threat  has  amplitude  fluctuations  (always  assumed  Rayleigh),  then  h  (f,t) 
is  a  zero  mean  nomally  distributed  random  variable  with  an  aatocavarianee 
defined  by 

h*(f.T)hs'(f,,,T,T  ■  i(f-f‘)«(T-r‘)f2(f.T)  (Sa) 

»  0  (5b) 

where  F,(f,t)  is  the  generalised  power  spectrum  and  l(x-x‘)  is  the  i>*rac 
delta  function  defined  for  any  function  f(x)  by 

d*  S(X-X*)f(x)  «  f(x') 

The  generalised  power  spectrum  is  parameter* ted  by  ^  ,  the  scintillated 
signal  decorrelation  time,  and  f  »  the  frequency  selective  bandwidth. 

There  are  two  forms  for  the  generalised  power  spectrum  depending  on  the 
product  f0T  where  T  is  the  minimum  sjrabet  period.  For  £  TM 


r2(f.r) 


at 


(6) 


1.864T  <S(t) 

0 

T  7  2*l~7 

1*S.S72t ‘f  “ 

l  o  J 


This  represents  the  flat  fade  condition  with  rospect  to  the  symbol  period, 

T.  for  f  T-l 
o 


r,(f,T)  •  25/4T»1/2^--expj-  — j|(iTfTfl)2.2itf’Tj2-(itfTo)2 


M  1 

C  1 

• 

Ids  exp 

J 

4.2 
-x  -2x 

«CP 

4  '  Cj  '  7/ 

where  P  =»  fa(l*C2)'J 

Cj  s  delay  parameter  (*0.25> 

for  both  equation  6  and  7 


(») 


Equation  7  provides  for  seint illation  random  delay  or,  equivalently, 
frequency  selective  effects. 

if  the  scintillation  threat  consist*  of  phase  effects  only  or 
if  there  are  no  scintillations  at  all,  thee 


S 


(9) 


*#(.\T)  -  5(0s(t) 

If  there  are  scintillations  of  any  kind,  then  there  5*  a  random 
time  varying  component  in  N(t)  in  addition  to  the  very  large  s  ■;  le 
(slow)  variations. 

The  total  electron  content  is 


N(t) 


Vt)j‘ 


df  g(f)exp(-i2itft) 


(Id) 


where  N{(t)  is  the  targe  scale  (slow)  component  and  g ( f )  is  a  tero  mean 
normally  distributed  variable  representing  the  random  component.  The 
autocovariance  of  g(f)  is 


where 


f<f 

*  *  *  *  t3 


eCf)  5  g*(’f) 


fg  *  Sayleigh  frequency 


(It) 


k'ith  the  specification  of  K  ■?’  h(t).  r  .  t'  f  .  and  ;(t), 

A  v  5  K  6 

we  obtain  a  complete  description  of  the  disturbed  siga»i  as  determined  fey 
the  propagation  environment  and  system  geometry.  £(b)  is  supplied  fey  the 
user  iadepeadeat  of  the  earireaseat  sad  geometry. 
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3. 


SIGNAL  SPECIFICATION  DEFINITION 


In  principle,  a  complete  signal  structure  specification  would 
include  all  of  the  time/ space  variations  possible  for  the  parameters  and 
functions  in  Equation  2.  In  practice,  it  is  sometimes  sufficient  to 
specify  extremum  values  or  simple  functions  for  the  propagation  quantities 
or  their  derivatives  accompanied  by  application  rules.  The  following  table 
lists  a  minimum  set  of  specification  parameters. 

TABLE  1 .  SIGNAL  SPECIFICATION  PARAMETERS 


{ 

f 

i 

i 

i 

h 


SPECIFICATION  PARAMETERS 


Maximum  Values 
Absorption,  KA 

2 

Energy  Angle  of  Arrival  Variance,  oQ 
Transmitter/Receiver  Vehicle 

dn 

Dynamics,  — —  ,  n=0,3 
dtn 

Signal  Decorrelation  Time,  tq 


Rayleigh  Frequency,  fD 

K 

Minimum  Values 

Signal  Decorrelation  Time,  tq 
Frequency  Selective  Bandwidth,  f 


10 


Two  nearly  universal  application  rules  apply  to  tq  and  £q.  For  r  , 

where  both  a  minimum  and  maximum  is  specified,  the  system  must  operate 

over  all  intermediate  values.  Similarly,  when  a  minimum  f  is 

specified,  the  system  must  handle  all  f  from  the  minimum  to  the  carrier 

frequency.  Exercising  these  ranges  are  necessary  because  tho  maximum 

performance  degradation  may  not  occur  at  tho  extreme  values  of  either 

t  or  f  . 

0  0 


Tho  parameters  to  be  specified  are  functions  of  tho  carrior 
frequency,  the  propagation  scenario,  the  link  gcoraetrios,  and  the 
velocities  of  the  system  segments.  If  possible,  the  specification 
should  cover  all  of  the  possible  signal  variations  realizod  by  exercising 
the  above  factors  ovor  their  entire  range.  Meeting  this  specification 
with  each  link  independently  would  provide  the  utlinmte  in  survivability, 
that  is,  a  system  that  can  survive  any  scenario  or  circumstance.  If  this 
specification  cannot  bo  met,  then  a  less  severe  specification  may  be 
possible  at  some  acceptable  loss  of  performance.  For  example,  in  nuclear 
environments,  signal  specifications  are  a  strong  function  of  the  maximum 
acceptable  outage  time.  Accepting  longer  outages  can  often  provide 
significant  specification  relief.  Additional  relief  might  be  possible 
for  systems  that  have  multiple  links  that  penetrate  the  ionosphere  at 
widely  dispersed  points.  Because  of  the  inherent  space  diversity, 
the  specification  would  not  have  to  cover  the  most  severe  case  for 
any  link,  but  some  reduced  level  of  threat.  Regardless  of  the  reason, 
however,  any  reduction  in  the  specification  results  in  some  loss  of 
system  applicability. 

Methods  to  calculate  the  propagation  parameters  are  described 
in  reference  1,  excluding  absorption,  which  is  covered  in  reference  2. 
Calculations  based  on  these  methods  serve  as  the  basis  for  developing 
signal  specifications.  The  first  step  in  the  specification  development 
process  is  to  choose  the  throat  environments.  For  natural  environments, 


the  threat  might  be  either  the  equatorial  or  polar  ionosphere  during  the 
solar  maximum  depending  on  the  location  of  the  system  links.  Nuclear 
threat  environments  would  be  calculated  by  computer  codes  specifically 
designed  for  that  purpose  from  plausible  burst  yield,  altitude,  and 
location  combinations.  The  next  stop  is  to  find  the  possible  link 
geometries  that  maximize  the  system  degradation  effects.  This  usually 
moans  minimizing  the  angle  between  the  link  lino  of  sight  and  the  earth's 
magnetic  field  where  the  environment  is  most  severe  and  maximizing  the 
link  path  length  through  the  disturbed  propagation  medium.  Those  two 
criteria  can  usually  be  satisfied  simultaneously  but,  if  not,  the  first 
almost  always  takes  precedence.  The  propagation  methods  in  references 
1  and  2  are  then  used  to  calculate  the  required  propagation  qua."* ties 
as  a  function  of  carrier  frequency,  area  coverage,  and  time.  The  system 
segment  velocities  are  chosen  to  provide  maximum  and  minimum  values  of 
the  signal  decorrelation  time.  Area  coverage  is  typically  represented 
in  contour  maps  projected  on  the  earth  of  the  propagation  quantities 
which  simultaneously  show  the  area  coverage  for  multiple  values  of  those 
quantities.  In  ambient  environments,  contour  map  sets  might  be  generated 
for  different  probabilities  of  occurrence  or  different  times.  The  nuclear 
environment  contour  sets  would  correspond  to  discrete  times  relative 
to  the  times  of  buvst. 

These  contours,  whether  ambient  or  nuclear,  provide  a  data  base 
for  the  final  specification  development.  The  preferred  specification 
would  refleet  the  most  severe  conditions.  If  this  is  not  feasible,  then 
the  propagation  data  would  permit  the  necessary  trade  offs  between  the 
threat  severity,  the  applicability  and  practicality  of  the  system. 

The  final  result  should  be  a  specification  that,  if  met,  would  provide 
adequate  system  performance  at  acceptable  cost  and  technical  risk. 


4. 


SIGNAL  SPECIFICATION  IMPLEMENTATION 


The  disturbed  signal  effects  represented  by  Equation  2  and 
the  specification  parameters  are  applied  twice  during  a  typical  system 
acquisition.  First,  the  signal  specification  is  used  during  system 
design  to  evaluate  candidate  design  solutions.  This  evaluation  is 
frequently  done  using  computer  digital  simulation  and  modeling  techniques 
which  provide  both  simulated  degraded  signals  and  detailed  descriptions 
of  the  dynamic  behavior  of  the  candidate  designs  over  all  the  specified 
signal  conditions  rapidly  and  economically. 

The  signal  specification  is  also  applied  during  testing  and 
evaluation.  At  a  minimum,  the  specification  determines  the  range  of 
effects  over  which  the  system  must  be  tested.  At  the  other  extreme, 
the  specification  may  define  not  only  the  parameter  ranges  but  also 
the  exact  form  of  the  signals.  This  latter  circumstance  occurs  when 
a  system,  in  its  full  operational  configuration,  cannot  be  exereised  in 
the  maximum  threat  environments.  For  example,  the  most  severe  natural 
scintillations  are  coincident  with  the  solar  sunspot  maximum  which 
occurs  only  once  every  eleven  years.  Systems  that  must  operate  in  nuclear 
degraded  environments  cannot  be  tested  at  all  because  of  the  1963 
Atmospheric  Test  "in  Treaty.  In  these  instances,  the  best  remaining 
method  is  to  degrade  the  satellite  signals  artificially  by  appropriately 
designed  link  simulators.  Ey  placing  a  simulator  at  one  or  more  points 
in  the  operation?!  system,  testing  can  be  done  limited  only  by  the 
accuracy  of  and  confidence  in  the  signal  specification. 

WTten  lint  simulators  are  necessary  for  testing  and  evaluation, 
questions  arise  over  what  constitutes  a  necessary  and  sufficient  test 
program.  For  example,  is  it  necessary  to  include  all  of  the  effects 
represented  by  Equation  2  in  every  test  .“'itsultaneously?  Is  it  necessary 
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to  test  simultaneously  every  link  in  a  system  or  even  every  link 
individually?  How  do  you  test  a  system  to  include  the  responses  of  the 
operators  as  they  react  to  degradation  on  system  links  from  propagation 
or  other  simultaneous  threats?  On  a  longer  term,  is  repeated  testing 
necessary  to  insure  that  the  system  hardware,  the  operational  procedures, 
and  the  operator  responses  remain  adequate  to  maintain  the  required 
performance?  Those  and  other  similar  issues  need  consideration  for 
each  system  tested.  The  answers  to  these  questions,  broadly  known  as 
"compliance  standards,"  should  accompany  the  signal  specification  for 
each  system  to  provide  a  completely  defined  test  and  evaluation  program. 

5 ,  SUMMARY 

The  preceding  sections  have  detailed  the  form  of  a  signal 

3 

structure  specification  for  application  to  satellite  C  systems.  The 
specification  is  intended  as  a  design  tool  as  well  as  a  definition  of 
the  signal  conditions  in  which  the  system  must  operate.  Adequate 
techniques  exist  to  apply  the  specification  both  in  computer  simulation 
during  design  and  by  3  link  simulator  for  test  and  evaluation. 

Appendix  A  describes  methods  to  implement  specific  realizations 
of  the  channel  impulse  response  function.  Appendix  t>  provides  a  brief 
description  of  link  simulators. 
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APPENDIX  A 

THE  GENERATION  OF  DIGITAL  REALIZATIONS 
OF  THE  CHANNEL  IMPULSE  RESPONSE  FUNCTION 


The  channel  impulse  response  function  is 


h(t,T)»A(t)exp 


T 


cr  N(t)  C 

i~ - J  d(Af)  h^t.Af) 

C  *• 


where 


A(t)  *  absorption  and  antenna  loss 

i  a  carrier  frequency 

c  «  light  speed 

z(t)  «  propagation  path  length 

r  *  classical  electron  radius 

o 

N(t)  a  N?  (t)*NR(t) 

Nj It)  =  slow  trend  total  electron  content 
NR(t)  *  random  total  electron  content 


The  generation  of  digital  real! tat  tons  of  Equation  A-l.  in 
general,  has  three  steps,  first.  A(t),  N^Ct).  and  tit)  are  discretized. 
Next,  samples  of  h  tt.Af)  and  Nett)  are  calculated  using  Monte  Carla 
techniques,  filially,  the  integral  over  Af  is  evaluated.  Most  of  the 
complexity  in  generating  these  digital  representations  lies  in  the  random 

t- 

sampling  of  h^U.Af)  and  N^U).  Most  of  this  appendix  will  concentrate 


on  these  sampling  techniques. 


*>/ 

The  statistical  properties  of  l»s(t,Af)  are  defined  by 


h  ( t ,  A  f ) 
s 


S 


te*  »o 


df  h  ( f,t )exp( i2«Aft*i2ii 


ft) 


(A- 2) 


where  h  (f,t)  is  a  zero  mean  normally  distributed  random  variable  with 
its  autocovariance  defined  by 


h*(f.t) 

h.(f\T»)  *  5  (f-f*)5(r-t*)  1'  (f,T) 

(A- 3a) 

hs(f.t) 

XCvJn  *  o 

(A- 3b) 

where  5(f)  »  Dirac  delta  function 

r,(f,t)  «  generalized  power  spectrum 

There  are  two  distinct  cases  to  be  considered,  first,  for 
flat  fading  defined  by  f  T*l  where  f  is  the  frequency  selective  band- 
width  and  T  is  the  symbol  period,  the  generalized  power  spectrum  is 

J.SMt  HO 

I\  ( f .  t )  5  -1  "  -=  wr  t A*  4 ) 

h*B.s?jtfowr 

where  ?  is  the  signal  decorrelation  time.  Because  of  the  delta 
correlation  in  deiay.  h^tt.if)  is  not  a  function  of  if  and  the  problem 
is  reduced  to  generating  a  sequence  that  ts  eaiy  a  function  of  time. 
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For  any  Lf,  let 


IUt.,Af)  *  p?i(t.  ^AtXi^V3-.  (A-5) 

SI  S  l  •  4  l 

>  t, 

u.  »  P  u.  J  ♦  u-p")  J  (A-6) 

where  p  a  cxp(-.M46  Al/x  ) 

o 

^“V'i-l 

g.  «  tevo  mean  normally  distributed  complex  sample 

gTgV  s  tanli  (2.146At/to) 


The  sequence  can  be  initiated  by  choosing  an  initial  value  tor  U.lt.,Af) 

*  * 

of  order  unity  and  a  value  for  u,  of  order  g* g. 1  .  Appendix  e  describes 
a  simple  algorUhm  for  sampling  zero  mean  normally  distributed  complex 
variables.  The  samples  from  fquation  A.$  represent  instantaneous  samples 
and  no  filtering  is  implied.  Also,  the  following  conditions  should  be 
met. 


f  /Atilo  (A«?a) 

M^lOd?  /At  (A* 7b) 

*  *b 

where  M  is  the  minimum  austber  of  points  in  any  sequence  generated  by 
negation 
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4N 


1 


H 


§m 


The  second  case,  defined  by  f  T<1,  is  frequency  selective  fading, 
liquation  A-2,  in  finite  difference  form,  is 


M 


hs(ti*AtVaAtV-MAt  ^  VVV°Xp(i2’Wi2’AVj)(A-8) 


where 


Af*  «  l/(2MAt) 

s>  ■  'r'i-t 

f  =  nA  f  * 
n 

t,  *  iAt,  | i | <M 
|Afj,|  <  1/(2At) 


From  Equations  A- 3a  and  A- 3b, 
be  written  as 


the  statistics  of  ^ )  can 


?  f.*At/2 


h*(f  ,t,)h  (f  ,t . )  @  8  <s. 

s  n  r  s '  a  i  tan  uIai 


t .-At/2 


/ 


f„*A  t*n 

df 

VAf,/2 


r,lf.O 


tA-@a) 


*  0  {A^b) 

^T?“Ttp  *  o  (A-9c) 
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I 


where,  for  f  T<1, 
0 


f  » 

Cj  a  delay  parameter  (”l).2J») 

Equation  A- 10  and  related  quantities  are  evaluated  in  Appendix  C. 
Generating  the  h  (t.,Af.)  consists  of  evaluating  the  h.(f  ,*,)  using 
liquations  A-9,  A-10,  and  Appendix  B  and  then  evaluating  Equation  A-8. 
for  adequate  statistical  sampling  and  numerical  resolution,  it  is 
necessary  that 


T  /At  5.  It) 

0 

(A- 1  la) 

K>5tH  /At 

(A- lib) 

A*«T/5 

(A.  Uc) 

A(Af)<l/aasU.l/f  ,U>t) 

Q 

(A-ild) 

[t  f|s0.7S 

6 

(A* lie) 

-0.2S<2»f  uS.45 

(A* Ilf ) 

where  (Af)»Af^-Af  .  Condition  A-llb  is  a  bare  minimum.  Larger  M  or 
multiple  sequences  are  advisable,  liquations  A-lle  and  A-llf  prescribe 
the  ranges  of  t  and  f  in  liquation  A- 10  that  include  95  percent  the 
total  energy.  The  coefficients,  h^tf^.tj),  are  not  necessary  outside 
those  ranges,  liquation  A-llg  and  the  T  dependent  tern  in  liquation  A-ild 
reflect  an  estimate  of  the  required  delay  range  to  handle  dispersion. 

Hie  integral  in  liquation  A-9  can  be  evaluated  to  one  percent  with 
trapezoidal  integration  with  frequency  and  delay  increments  less  than  or 
equal  to  0.03/t^  and  0.03/f  ,  respectively. 

The  representation  used  for  the  random  portion  of  the  totai 
electron  content  depends  on  the  application.  The  first  use  of  N{  ( t )  in 
liquation  A-l  is  in  the  carrier  phase  where  it  can  degrade  phase  or 
frequency  acquisition  and  tracking,  for  plus©  effects 


M* 

Mf  .IsAf"  £  g(f’)exp(.i2itf*t.  l  (A'  12) 

n=-M*  1  1 

where  Af“  =*  T / 1 -M  *  a  t  * ) 

£'’  *  'r'i-t 

f  *  =  nt  f  “ 

a 

1 1  *  Ut\  j  i|fM; 

g(f^)  is  a  sere  normally  distributed  variable  whose  resaisiag 

statistics  are  defined  by 


g-(fj)gTf;> 


f  .f ' 
;v  I** 

’)  4*-  Af"/2 


where  g(f*  )  a  »*(**  } 

*  -n  K  '  n 

2  /roc't.m\ 

*  \T'/ 

fR  *  Rayleigh  frequency 

initial  estimates  for  At'  and  M*  are 

A«*<i/(4fR) 

M*»SU 

for  frequency  tracking.  Equation  A-U  cannot  he  used,  instead,  let 


TA(fc/roc)2 

p^7T?p  ■ 

(A-13a) 

-  0 

K;K 

(A-I3b) 

(A.  Ha) 
CA-i4b) 


VV  *  *,:„."**  8^a>tP*pU3»^4* 


(A«i5) 


mi mm* 


The  second  use  of  N(t)  in  Equation  A-l  is  in  the  dispersive 
phase.  NR(t)  can  usually  be  omitted  because  the  dispersive  effects  are 
primarily  determined  by  the  magnitude  of  N(t)  and  not  the  time  dependence. 

The  last  use  of  NR(t)  determines  the  group  delay  and  group  delay 
rate  which  can  degrade  time  synchronization  if  sufficiently  large. 
Equations  A- 15  and  A- 16  should  be  used  with  fR  replaced  by  f^.. 


i 

| 

i 

] 

a 

I 

I 


W  [■‘(lM'VJ)2]  h 


(A- 17)  | 


If  simultaneous  handling  of  the  total  electron  content  effects  is 
necessary,  then  Equations  A- 15  and  A- 16  ai'e  necessary.  The  substitution 
of  should  not  be  used.  Unfortunately,  the  number  of  time  increments 
necessary  to  span  the  time  range  can  be  very  large.  Exercising  the  range 
of  time  may  not  be  necessary  in  all  cases.  Some  experimentation  is  often 
indicated.  Another  problem  with  Equation  A-15  is  that  the  sura  cannot  be 
evaluated  with  fast  fourier  transforms.  Finally,  multiple  sequences  of 
N^t^)  should  be  used  to  test  the  adequacy  of  the  sampling  for  either 
Equation  A- 12  or  A-15. 


The  entire  channel  impulse  response  function  is 


2rf  2(t  )  cr  N(t  )  l  % 

h(t.  ,t  .)=A(t.)exp  i - - — i--i - j -  A(Af)  E  hs(t.,Afj.) 

c  k= **  L 


crN(t.)Af^  .  /  Z<V  croN(V 

exp  .i - -  l2irAfk(  y  c  -  2v{2 


(A-18) 


m 


where 


1=1/ (2AtA (Af) ) 


cr  N(t , ) 
01 


<Lflr 


In  principle,  one  only  needs  to  insert  the  proper  quantities 
discretized  to  the  smallest  increment  in  time,  delay,  or  Af^  and  sum 
over  Afj..  In  practice,  however,  the  wide  range  of  the  delay  and  time 
necessary  to  adequately  sample  the  various  functions  make  this  impractical. 

It  is  usually  necessary  to  handle  some  of  the  terns  independently, 
particularly  the  total  electron  content  and  path  length  determined  effects. 

The  channel  impulse  function  is  thereby  greatly  simplified. 

h(t.,y»A(t.)hs(ti,r.)  IA-19) 

where  for  f  Tsl 

O 

M 

h  (t . ,t . i®Af '  £  h  {f  t  }exp(i2*f  t.) 

*  1  J  n@-M  *  *  *  ”  * 

or  for  f  T> 1 

(5 

hs(ti‘V*VVAf)*  Tj8° 

«U  ,  ijjo 

A  fortran  program  for  generating  h  (t.,?.)  can  be  found  in  Appendix  E. 

Equation  A- 19  is  impleaentable  in  a  link  simulator  as  described  in  Appendix  0. 

An  important  consideration  for  implementing  Equation  A-J8  or  A- 19 
is  the  dynamic  range  required  to  resolve  the  amplitude  fluctuations,  if 
present.  The  amplitude  of  the  channel  impulse  response  function  at  each 
discrete  delay  is  a  Rayleigh  distributed  variable.  Let  the  amplitude  of 
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h(tj.Tj)  be  R.j.  The  distribution  function  for  R. ,  is 


2R.  .  /  f  o\ 

P(R_)a_U  oxp^-RTj/o'J 


(A-20) 


where 


t .♦At/2 


liquation  A-20  is  easily  invertable  to  form  probability  statements  about 


.  Thus, 


l -exp 


■(fl 


» probability  that  R.^<Rr 


"  probability  that 


We  require  that  the  amplitude  not  exceed  the  dynamic  range  99.8  percent 
ot  the  time.  The  upper  and  lower  bounds  on  the  amplitude  are, 
respectively 


li/Vu  5  5  2.6S 

3  £*in(0.999)j  2  *  5.16; 


16x10  * 
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where  the  probabil ity  that  the  amplitude  exceeds  the  dynamic  range,  0.002, 
is  equally  split  between  the  limits.  In  decibels 


3OdB<20 


Thus.  38.4  dB  of  dynamic  range  is  required.  If  this  is  not  feasible,  then 
the  loss  of  accuracy  should  be  at  the  largest  amplitudes.  This  is 
because  system  performance  is  much  move  sensitive  to  amplitude  fades  than 
to  amplitude  enhancements . 

The  resolution  in  the  amplitude  and  phase  is  a  function  of  the 
resolution  of  the  real  and  imaginary  components  of  h(t^,tj).  The  poorest 
resolution  occurs  in  deep  fades  at  the  smallest  amplitude.  The  maximum 
increment  is  defined  by  requiring  that  successive  values  of  either 
component  of  hlt.,tj)  not  vary  by  more  than  ten  percent.  This  gives 
adequate  accuracy  on  the  amplitude  and  resolves  the  phase  to  within  six 
degrees,  tf  the  component  increment  is  fixed,  then  the  smallest  value  must 
be  used.  The  fixed  increment  for  either  component  can  be  expressed  as 
one  tenth  of  (R^/e^.  Variable  increments  are  more  desirable 
because  they  permit  more  efficient  storage  of  component  data.  They  also 
permit  the  use  of  logarithmic  digital  to  analog  converters  in  link 
simulators  as  described  in  Appendix  b. 

The  total  electron  content  effects  neglected  in  liquation  A*  19 
can  be  bounded  to  determine  necessary  but  not  sufficient  conditions  for 
successful  systt-«  operation.  The  total  electron  content  and  its 
derivatives  are  gaussian  distributed  random  variables  with  the  following 
means  and  standard  deviations. 
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■-.n*?  «TS^8sti 


dmN(t) 

dtm 


a 


(A-21a) 


o 


No  3 


(2.f)2,,o(fc/r0c)2 

‘p.(2„>7]S2 


(A-21b) 


These  statistical  moments,  when  converted  to  moments  of  phase,  doppler, 
dopplev  rate,  jerk,  etc.,  allow  performance  degradation  estimates.  These 
statistics  also  determine  the  amount  of  time  that  any  of  the  derivatives 
exceed  some  value.  As  such,  they  provide  the  means  for  specifying  various 
system  parameters  consistent  with  some  performance  requirement. 


The  total  electron  content  drives  the  dispersive  distortion  of 
the  signal  waveform.  As  a  measure  of  this  distortion,  we  use  the  expec¬ 
tation  value  of  the  absolute  value  of  the  delay  coordinate. 
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er  N(t)Af* 

I  dt  [  r  f 

d(af)exp 

-  -i2»Aft 
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ftp 

c  J 

(A-22a) 
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er  N(t) 
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Since  the  statistical  moments  of  the  measure  of  distortion  are  proportional 
to  the  moments  of  the  total  electron  content,  performance  criteria  can 
be  developed  as  a  function  of  modulation  and  other  system  details.  The 
criteria  should  guarantee  that  the  ratio  of  ;?|  to  T  be  small. 


The  finite  difference  equations  have  been  written  in  the  lowest 
order.  This  was  done  for  simplicity  and  to  be  consistent  with  fast 
fourier  transform  algorithm-. .  These  algorithms  can  evaluate  liquations 
A-8,  A-12,  A-18,  and  A-19  much  wore  efficiently  than  conventional  methods. 
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APPENDIX  B 

A  RANDOM  NUMBER  GENERATOR  FOR 
NORMALLY  DISTRIBUTED  COMPLEX  NUMBERS 


Let  g  be  a  normally  distributed  complex  number  defined  by 


-  2 

8  «  (tl-la) 

S8  B  0  (B-ib) 

S  *  «  (B-Ie) 

P(|«|)  •  -jr W(’UI 'fa")  (B-ld) 

«• 

0 


where  P(lgj)  is  the  probability  distribution  for  the  amplitude  of  g.  Let 
be  a  uniformly  distributed  random  variable  on  the  interval,  ixR^l. 

A  sample  of  g  is  calculated  by 


S  9 


18-2) 


APPENDIX  C 

EVALUATION  OF  THE  GENERALIZED,  DELAY. 
AND  FREQUENCY  POWER  SPECTRUHS  FOR  f  T<1 


The  generalized  power  spectrum  is 


where  f  *  f.UKp'4 

(.‘j  5  delay  parameter  (=o.JS) 
f0  =  fluency  selective  bandwidth 
fQ  5  sisaal  decarrelat ion  time 

A  m  te  ^uatiea  C-l.  geed  te  an*  percent,  is 
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i(r^) 


t  ?  {i*O.4?04?f 


a  fff'expi 


(2lff'T) 
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2C 


*  Ct )  .  2*0 


(C-Sa) 


r,(t)  «  «f'exp|*-  -2irf*t^|2-P(t)exp(-22)J  ,  2^0 

Ht)  a  ^(0.7478SS6t -0.0058798)  t  ♦  0.5480242^6 

liquation  C-5  is  accurate  to  about  five  figures.  The  frequency  power 
spectrum  is  defined  as 


(C-Sb) 


u* 

r2(f)  * Jdtr^if.t) 

t?(0  *  »}tQ«p^(»tQfj*  J 


(C-h) 


(C-7) 


The  normalisation  integrals  for  the  above  spectrum*  are 


«-  a. 


H 
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df  -  l 


(C»8) 


it)  »  1 


(C-9) 


dfr\tf)  *  l 


ss 


Frequently,  integrals  over  the  various  power  spec t rums  are 
necessary  as,  for  example,  in  liquation  A -9.  Trapezoidal  integration  with 
the  following  integration  increments  is  accurate  to  approximately  one 
percent . 

afj  i  0.03/to  (C-lla) 

Atj  s  O.O3/f0  (C-llh) 


APPENDIX  D 

PROPAGATION  EFFECTS  LINK  SIMULATORS 


Figure  D- 1  shows  a  general  implementation  of  a  link  simulator 
which  is  the  hardware  analog  of  the  channel  impulse  response  function 
defined  by 


U(t)  «|  dthft.t )S(C-t) 


(P-l) 


II  ! 


wnere  J>(t)  -  transmitted  (input)  signal 

R(t)  «  received  (output)  signal 
htt.T)  *  channel,  impulse  response  function 


the  link  simulator  represents  a  discretized  version  of  R-l  <4*e?e 


yti  * 


BMjaealjhtyy^ 


(&*:a) 


Q^Ct)  *  R/A|faag|hCt ]| 


(&-2M 


R/a{  }  3  digital  to  analog  conversion 

Algorithms  to  calculate  values  of  tne  ekas**!  impulse  response  fuactian 
were  detailed  ia  Appendi*  A.  the  sampling  respai reseats  in  Append is  A 


F.gure  D-l .  Propagation  Effects  Simulator 


are  sufficient  to  support  the  digital  to  analog  operation  without  any 
additional  filtering  to  suppress  aliasing. 

Figure D-2  shows  a  particularly  simple  link  simulator  for  flat 
fading  and  when  all  other  effects  are  ignored  or  simulated  separately. 

The  delay  is  used  to  get  an  independent  noise  source  for  the  I 
channel  without  an  extra  noise  generator.  The  filters  are  single  pole 
providing  in  cascade  the  power  spectrum  represented  by  Equation  6.  For 
a  RC  filter,  the  time  constant  is 


RC  =  t  /2.146  CD-3) 

o 


The  input/output  can  be  at  any  frequency  level. 


APPENDIX  E 

CIRF,  A  FORTRAN  PROGRAM  FOR  GENERATING  SAMPLE  DIGITAL 
REPRESENTATIONS  OF  THE  CHANNEL  IMPULSE  RESPONSE  FUNCTION 


This  appendix  contains  a  fortran  program  to  sample  digital 
representations  of  the  channel  impulse  response  function  in  Equation  A-19 
less  the  absorption  and  antenna  loss  terms.  The  coding  is  consistent  with 
Fortran  IV  with  the  exception  of  the  dimension  declaration  in  the  fast 
fourier  transform  subroutines  (see  comment  in  SUBROUTINE  SETD  list). 

The  methods  used  to  generate  the  frequency  selective  sequences 
are  also  used  for  the  flat  fade  sequences  rather  than  Equations  A-5  and 
A-6 .  The  flat  fade  sequence  generation  is  thus  slower,  but  the  sequences 
generated  are  continuous  from  the  last  sequence  term  back  around  to  the 
first  term,  a  very  useful  property  for  hardware  testing  applications. 

The  flat  fade  sequences  are  nevertheless  consistent  with  Equation  A-4. 

The  program  tests  the  input  for  compliance  with  Equations  A-lla, 
A-llb,  A-lle,  and  A- Ilf.  The  recommended  integration  increments  are  also 
enforced.  The  use  of  floating  point  variables  for  the  channel  impulse 
response  function  output  guarantees  compliance  with  the  stated  phase  and 
amplitude  resolution  requirements  on  all  known  computer  systems  that 
support  Fortran  IV.  The  user  is  responsible  for  satisfying  Equation  A-Uc. 

The  code  provides  time  sequences  of  the  impulse  response  function. 
The  first  sequence  centers  on  a  delay  of  -O.04/F0+PD/2  where  FO  is  the 
frequency  selective  bandwidth  and  DD  is  the  delay  increment.  The  main 
program  loop  over  delay  ends  at  statement  label  70.  Just  prior  to  that 
statement,  the  time  sequence  for  the  current  delay  is  in  the  H  array. 
H(I*1-1)  and  H(I*1)  represent  the  in-phase  (real)  and  quadrature  (imaginary) 
impulse  response  function  components,  respectively.  The  user  is 


responsible  for  providing  the  extra  coding  for  disposing  of  the  sequence 
data.  The  data  must  be  moved  out  of  the  array  before  the  loop  returns  to 
calculate  the  next  sequence. 

Any  remaining  questions  should  be  referred  to 
Director 

Defense  Nuclear  Agency 
RAAE/Major  Leon  A.  Wittwer 
Washington,  DC  20305 

Autovon  221-7028 
Commercial  202/325-7028 

The  program  was  programmed  and  debugged  on  a  Heathkit  H89 
Computer  system. 


PROGRAM  CIRF 
C 

C  THIS  PROGRAM  CALCULATES  RANDOM  SAMPLES  OF  THE  SCINTILLATION  CHANNEL 
C  IMPULSE  RESPONSE  FUNCTION  AS  DEFINED  IN  "A  TRANS-IONOSPHERIC  SIGNAL 
C  SPECIFICATION  FOR  SATELLITE  C**3  APPLICATIONS* “  (ENCLOSING  DOCUMENT) 

C 

C  THE  PROGRAM  INPUTS  ARE 
C  TAUO*  SIGNAL  DECORRELATION  TIME (SEC) 

C  FO,  FREQUENCY  SELECTIVE  BANDWIDTH (HZ) 

C  DT,  TIME  INCREMENT (SEC) 

C  (THE  TIME  INCREMENT  MUST  BE  LESS  THAN  OR  EQUAL  TO  TAUO  DIVIDED 

C  BY  TEN) 

C  DD.  DELAY  INCREMENT (SEC) 

C  (THE  DELAY  INCREMENT  SHOULD  BE  LESS  THAN  OF  EQUAL  TO  THE 

C  RELEVANT  SYSTEM  SYMBOL  DIVIDED  BY  FIVE.  IF  DD  IS  GREATER  THAN 

C  OR  EQUAL  TO  SIX  TENTHS  DIVIDED  BY  FO.  THEN  ND ( SEE  BELOW)  IS 

C  SET  TO  ONE  AND  THE  FLAT  FADE  SPECTRUM  IS  USED) 

C  NT,  NUMBER  OF  DISCRETE  TIMES 

C  (NT  MUST  BE  A  POWER  OF  TWO.  NT  MUST  BE  SUFFICIENTLY  LARGE  SUCH 

C  THAT  NT  TIMES  DT  IS  GREATER  THAN  OR  EQUAL  TO  ONE  HUNDRED  TIMES 

C  TAUO) 

C  NO.  NUMBER  OF  DISCRETE  DELAYS 

C  (ND  MUST  BE  SUFFICIENTLY  LARGE  SUCH  THAT  ND  TIMES  OD  IS  GREATER 

C  THAN  OR  EQUAL  TO  SIX  TENTHS  DIVIDED  BY  FO) 

C  SEED.  SEED  FOR  RANDOM  NUMBER  GENERATOR 
C  (SEED  MUST  BE  AN  ODD  REAL  NUMBER) 

C  REQUIRED  FUNCTIONS/SUBROUTINES 
C  FLOAT  (FIXED  TO  FLOATING  POINT  CONVERSION) 

C  SETD  ( INITIALIZES  FAST  FOURIER  TRANSFORM  TABLE) 

C  IFIX  (FLOATING  TO  FIXED  POINT  CONVERSION) 

C  DINT  (DELAY  INTEGRATOR) 

C  OF AM  (PF*( STANDARD  DEVIATION  OF  FOURIER  COEFFICIENT)) 

C  RANCOF  (RANDOM  OF* (FOURIER  COEFFIC l ENT )  GENERATOR) 

C  PFT  AND  FDD  (FAST  FOURIER  TRANSFORM) 

C  POWER  (SIONAL  POWER  DENSITY  AS  A  FUNCTION  OF  DELAY) 

C  TCQR  (COMPLEX  SIGNAL  TIME  DECORRELAT ION  FUNCTION) 

C 

DOUBLE  PRECISION  SEED 

C 

C  THE  FOLLOWING  DIMENSION  LIMITS  ARE  MINIMUM  VALUES 
c  DIMENSION  H(2*Nf >.0(NT) 

C 

DIKE  NS  ION  H  ( 40*3*  > .  Pt  20*0  > 

COMMON  /DAT/NFP.PFP.NDP.&DP.ND.&O 
COMMON  /FSD/Ct . T AUO . FO 
COMMON  /RAfl/SEED 
COMMON  /BFM/BFMS 
C 

C  SET  DELAY  PARAMETER 

c 

DATA  eu.24^ 

e 


C  READ  AND  TEST  DATA 
C 

C  INPUT  IS  ON  LOGICAL  UNIT  6 
C  OUTPUT  IS  ON  LOGICAL  UNIT  2 
C  OUTPUT! CRT)  IS  ON  LOGICAL  UNIT  1 
C 

CALL  0PEN<6.'SYU INPUT.DAT  ') 

CALL  OPEN! 2* 'SYlt OUTPUT.  DAT  ') 

READ<6,1000>TAUO.FO,OT,DD,NT,ND.S£ED 

1000  FORMAT !4£10.2<2I5>D10.2) 

WRI TE ! 2 , 1001 >  T AUO . FO . OT .  DD . NT ,  NO , SEED 

1001  FORMAT! IX. IP, 27HS10NAL  DECORRELATION  TIME-  .ElO.2/  1X,31MFREQUENCY 
1  SELECTIVE  BANDWIDTH-  .E10.2/1X. 16HTIME  INCREMENT-  .E10.2/1X, 17HDE 
1LAY  INCREMENT-  .E10.2/1X . 17HNUMBER  OF  TIMES-  > 15/ IX, 16HNUMBER  OF  D 
l ELAYS-  .15/1X.6MSEED-  ,020.12) 

WRITE! 1,1010) 

1010  FORMAT ! 1 X «  3HRUNN I NO  ) 

IF! DT.LE. TAUO/ 10. )00  TO  10 
WR I TE (2> 1 002 ) 

1002  FORMAT 11 X.33MERR0R-TIME  INCREMENT  IS  TOO  LARGE  ) 

STOP 

10  AMP-FLOAT! NT) 

NP-0 

20  AMP-AMP/2. 

IF!AMP.OT. 1. >00  TO  20 
IFtAMP.EO. 1. )G0  10  30 
WRITE!?. 1003) 

1003  FORMAT! IX. 30HERR0R-NT  IS  NOT  A  POWER  OF  TWO  > 

STOP 

30  IF!FLOAT(NT)*OT.OE.100.*TAUO)00  TO  40 

WRITE!?, 1004) 

1004  FORMAT! IX. 35HERR0R- TIME  RANGE  LESS  THAN  100  TAUO  > 

STOP 

40  IFIDD*FL0AT!N0).0E..G/F0)00  TO  50 

WRITE!?. 1005) 

1005  FORMAT! IX. 34HEftR0R-D6LAV  RANGE  LESS  THAN  0.6/F0  ) 

STOP 

50  IF  t  DO.  LT , .  6/F01CQ  TO  60 

NO~t 

WRITE!?.  1004.) 

1006  FORMAT! IX. 39NFLAT  FADE  LIMIT-ONLY  ONE  DELAY  REGUIREO  ) 

60  CONT INUE 

C 

C  INITIALIZE  FAST  FOURIER  TRANSFORM  TADLE 
C 

CALL  SETDtO.NP) 

L*NT*NT-t 

C 

C  CALCULATE  FREOUENCV  INCREMENT  AND  INTEGRATION  VARIABLES 
C 

C  DRsFREOUENCY  INCREMENT 
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DF-1. /< FLOAT < NT )*0T> 

C  NFP«NUM6ER  OF  FREQUENCY  SUP INCREMENTS 
NFP*lFIXCJ3.<*TAUO*DF>-M 
C  PFP-FREOUENCY  SUP INCREMENT 
DFP-DF /FLOAT <NFP> 

C  NOP-NUMEER  OF  DELAY  SUB INCREMENTS 
NDP*IFIX(33.*DO*FO>M 
C  OOP-DELAV  SUP INCREMENT 
OOP-DD/FLOATINDO) 

c 

C  CALCULATE  TIME  SEQUENCES  FOR  EACH  DELAY 
C 

c 

C  DSTART  IS  THE  BEOINNINO  OF  THE  DELAY  WINDOW.  TWO  PER  CENT  OF  THE 
C  SIGNAL  ENERY  ARRIVES  BEFORE  DSTART  AND  THUS  IS  NEGLECTED  WHEN  DSTART 
C  IS  INITIALIZED  TO  -.04/F0. 

C 

DSTART*-. 04/F0 
C 

C  INITIALIZE  VARIABLES  FOR  STATISTICAL  TESTING.  THE  AVERAGE  SIGNAL 
C  POWER  AT  EACH  DELAY (THE  AVERAGE  VARIANCE  OF  THE  FOURIER  COEFFICIENTS 
C  TIMES  PD**2>  AND  THE  TIME  DECORRELATION  PROPERTIES  OVER  ALL  SEQUENCES 
C  NEAR  TAUO  WILL  BE  TESTEO. 

C 

I COR* IF I X  t  TAUO/OT ) 

!COR»ICOR*;CQR*l 

CORRsO. 

COR 1*0. 

DO  70  1*1. ND 
T*OSTART*.?.*D£> 

WRITE  1 1  *  lOli'IT 

1012  FORMAT! IT. IP. JtMGCNERATlNO  SEQUENCE  FOR  DELAY*  .EI0.2) 

W*L 

F?TART*.S*OF 
AMP* I. 

C 

C  ZERO  DC  FREQUENCY  COMPONENT  TO  MAINTAIN  ZERO  MEAN  SEQUENCE 
C 

HI  2>«Q. 

C 

c  DFMS*INITIAL  DELAY  INTEGRATION  POINT  FOR  DRAM 
C 

DFM?=.  S*OINT  f  fiSTART  .FSTART I 

po  eo  ytes.NT.s 

IFtAMP.LE.O.  >0©  TO  S'?. 

AMPsDFAMtDSTART.FSTAST f 
CALL 

CALL  P.ANCOF (Wt*  !.H(*,«ti.AMPl 
FSTAfiT?F-:TAftT*&F 

sa  to  eo 

OS 


<5 


HUM  1*0. 

H(K)»0. 

H<K*l )*0. 

80  K-K-2 

AHP-DFAM  <  PST ART .  FSTART ) 

CALL  RANCOFIHIKl.MIKMl.AMP) 

c 

C  use  FAST  FOURIER  TRANSFORM  TO  GENERATE  FINAL  SEQUENCE 
C 

CALL  PFT < H> D> 

CALL  PDD(H) 

C 

C  TEST  SIGNAL  POWER  IN  SEQUENCE  TERMS  AND  ACCUMULATE  DATA  FOR  TIME 
C  CORRELATION  TEST.  USE  TRAPEZOIDAL  INTEGRATION. 

C 

TVAR*0. 

DO  90  UMC0R.L.2 
K-U-ICORM 

CORR"CORR*H(K)*HU>*N(KM>»H(JM> 

CORI*CORl*M(KM 1  *MU1-H<K1*MUM ) 

90  TVAR=TVAf<*M  <  J  >  *M ( U ) *H I J* I >  *H<  J* 11 

AMP«2.0*TVAR*D0«»2/FL0AT<L-iC0R*2> 

PSUM* . 5*P0UER <  DSTAftT > 

T-DSTART 
DO  100  UM.NOP 
T®T*DDP 
SM«POWE»«T) 

100  PSUMbPSUM*SM 

PSUM* I  PSUM- . 5*$M ) *DDP 

T*D$TART*.S*PD 

UR  I  Tt:  ( 2 .  1 000 1  T .  AMP .  PSUM 

100C  FORMAT ( / 1  K .  IP. ?MOEl-AY=  .E 10.2/1 X . 24MAVSRA0F  SEGUCNC'E  POWER*  . 

IE10.2/U*  2JNTHE0RETICAL  SEQUENCE  POWER*  .£10.21 
C 

C  THE  H  ARRAY  NOU  CONTAINS  A  SAMPLE  DIGITAL  REPRESENTATION  OF  THE 
C  CHANNEL  IMPULSE  RESPONSE  FUNCTION  FOR  A  DELAY  OF  PSTAftT*. 5*PP.  THE 
C  IN-PHASE  AND  QUADRATURE  WJt.T {PLIERS (THg  REAL  AND  IMAGINARY  PARTS  OF 
C  THE  CHANNEL  IMPULSE  RESPONSE  FUNCTION.  RESPECTIVELY!  ARE  DEFINED  DV 

e 

C  IN=PHASE  COMPONENT  MULTIPLIER  •  MI  1*1-11 

e  QUADRATURE  COMPONENT  MULTIPLIER  *  M!I*t> 

e 

e  WHERE  I  IN.CXF.S  SUCCESSIVE  VALUES  IN  TIME.  THE  USER  MUST  DECIDE  ON 
C  THE  DISPOSITION  OF  Tm?  DATA  IN  THE  A  ARRAY  DEPOSE  PROCEEDING  ON  TO 
C  the  NE*T  DELAY.  USUALLY.  THE  DATA  IS  STOREO  AND  l.Af§R  USED  {N  AN 

C  integral  over  delay,  jn  these  integrations*  PD  must  t*£  USED 
o  explicitly  as  the  integration  increment. 
c 

>0  OSTAftT«OSTfiftT*PP 

e 

e  TEST  TIME  CORRELATION  PROPERTIES 

e 
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o  in  n  n  c>  n 


CORR-2.  0-CORR-DD--2/FLOAT  ( L  - ICOR-2  > 
CORI-2.0*C0RI*DD**2/FL0AT<L-IC0«»2> 

T-.S»FLOAT(tCGR-l)«DT 

AnP-TCOR(T) 

WRITE <2. 1  OOP  >  T  >  C  ORR i COR I • AMP 

1009  FORMAT (  / 1 X .  IP .  26HT1ME  DISPLACEMENT  TESTED-  .E10.2/  1X.26MCAI.CU4, ATE 
ID  DECORRELATION-  . 2E1 0.2/1 X. 25MThEORETICAL  CALCULATION-  .EI0.21 
WRITE (2. 1007) 

1007  FORMAT ( t  X • 10HEND  OF  RUN  ) 

END 

FUNCTION  OF AM (DSTART • FSTART > 

C 

C  THIS  FUNCTION  CALCULATES  DF  TIMES  THE  STANDARD  DEVIATION  OF  TME  RANDOM 
C  FOURIER  CORFF ICIENT  CORRESPOND I NO  TO  A  DELAY  OF  DSTART*. 3»DD  AND  A 
C  FREQUENCY  OF  FSTART*. S*DF.  THE  VARIANCE  OF  THIS  FOURIER  COEFFICIENT 
C  IS  THE  GENERALIZED  POWER  SPECTRAL  DENSITY  INTEGRATED  OVER  DELAY  FR.jM 
C  DSTART  TO  DSTART *DD  AND  OVER  FREQUENCY  FROM  FSTART  TO  FSTART*DF. 

C  INTEGRATION  IS  BY  THE  TRAPEZOIDAL  RULE. 

C 

C  ARGUMENTS 

C  DSTART- INITIAL  DELAY! SEC) 

C  FSTART- INITIAL  FREQUENCY (HZ ) 

C  REQUIRED  COMMON  VARIABLES 
C  NFP-NUMBER  C'F  FREQUENCY  SUB  INCREMENTS 
C  DFP-FREQUENC V  SUB I NCREMENT ( M2 > 

C  DDP-DELAY  SUB INCREMENT (SEC) 

C  DD-DELAY  INCREMENT 
C  DFMS-O. S*BINT ( DSTART . FSTART  > 

C  FUNCTIONS  REQUIRED 
C  DINT  (DELAY  INTEGRATOR) 

c  sort  (floating  point  square  rood 

c 

COMMON  /DAT/NFF‘.DFP.NtiP.DCiP.ND.DO 

COMMON  /UFM/DFM'J 

F-FSTART 

t»FAM=DFM'3 

DO  10  Jsl.NFP 

F«F*DFP 

SM^DINT ( DSTART. P) 

»0  pFAM?t<F6M*S« 

DFMt=.S*SM 

DF6M=S0.RT  <  &FP-DDP*(DF£M-OFMS » >/©€» 

RETURN 

END 

FUNCTION  DINT i ©S TART. P> 


THIS  FUNCTION  INTEGRATE?  THE  OF NEPAL. { JED  PQUER  SPECTRUM  FROM  &  ©ELLY 
OF  ©STAR?  TQ  ©START.©©  AT  A  FREQUENCY  ©E  F  AN©  ©{VIDCS  fttL  RESULT  ©Y 

BD*. 


ARGUMENTS 

©START-INITIAL  delay t SEC > 


uuu 


C  F-FREOU€NCY(MZ> 

C  REQUIRED  COMMON  VARIABLES 
C  NO-FLAT  FADE  FLAG 
C  NO- I. FLAT  FADING 

C  ND.OT. t. FREQUENCY  SELECTIVE  FADING 
C  NOP -NUMBER  OF  DELAY  SUB INCREMENTS 
C  DDP-OELAY  SUBINCREMENT 
c  FUNCTIONS  required 
C  PSDJ  {GENERAL  1 2EO  POWER  SPECTRUM) 

C  PSOFF  IFLAT  FACE  SPECTRUM) 

c 

COMMON  /DAT/NFP.OFP.NOP.DOP.NO.DO 
TEST  FOR  FLAT  FADING 
1F<NB.ECM)Q0  to  too 

0-PSTART 

DINT-.S*PS03(F.D> 

DO  10  I-t.NDP 
D-D-POP 
SM-PSP3 1 F .  0 1 
10  DINT-OINT-SM 

DINT-DINT-. 5*S« 

RETURN 

too  DINT-PSOFFtFI/OCR 
RETURN 
END 

FUNCTION  PSPT(F.D) 

C 

C  THIS  FUNCTION  is  THE  GENERALIZED  POWER  SPECTRUM  FOR  FREQUENCY 
C  SELECTIVE  CASES. 

C 

C  ARGUMENTS 
Z  F*FREOUgNCYCHJ> 

C  0=E>ftAV 

C  REWIRED  COMMON  VARIABLES 
C  Ct -DELAY  PARAMETER 
C  YAUO>$tONAL  DECORRELATION  TIME 
C  FQ-FREOUENCV  SELECTIVE  SA»£>WriOT* 
e  functions  required 
e  sort  {floating  point  square  root) 

c  AOS  IfLOaTSfiG  POINT  ABSOLUTE  VALUS  * 

C  EXP  (FLOATING  POINT  NATURAL  EXPONENTIAL) 

c 

COMMON  /PSSt/Ct.TAUO.FO 
FPRIH£=rFO*ioftT(  I.*Ct*#^) 

TPFPD-X-.  5*FPR  IME*D 

TO? i  *  *  c  t  »  H  S.  t  *s* T  AL^*F  >  **i‘- TPFRO » /C  t » 
|F(ADS(Z).LT.t.)GO  TO  40 

IFIf.LT.O.  )6-0  T6  JO 

o?*t .  »s*eip»-z**L  »*♦  t .  *  ( « .  t*T-.  t  t;*»x  )**-.  t?t  t*»  i/sqnt  t  zt 
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d  o  d  f »  n  o  n 


00  TO  30 

20  OZ-l.775*<l.*U.OO»4#X-.l?9)*X*.29>*X>/SaRT<-Z> 

00  TO  30 

10  02-( <  < .03»*Z-. J09)#Z-.729>*Z-.675>*2 

OZ-1.3I3*EXP<OZ> 

50  P$03»2.98l»FPfclME*TAu0«exP(.3*Cl**2-TPFPD»*GZ/SGftT«,i> 

RETURN 

END 

FUNCTION  PSDFF(F) 

C 

c  this  function  is  the  flat  fade  power  spectrum,  this  function  IS  NOT 

C  THE  GENERALIZED  POWER  SPEC TRUMt FUNCTION  PSD3>  INTEGRATED  OVER  ALL 
C  DELAV.  A  CLOSE  APPROX IHATtON  NAS  SEEN  CHOSEN  INSTEAD  BECAUSE  OF  ITS 

e  utility  in  software  and  hardware  applications,  white  noise  filtered 
C  BY  two  single  pole  filters  reproduces  this  spectrum. 
c 

C  ARGUMENTS 
C  F-FREOUENCYfMZ) 

C  REQUIRED  COMMON  VARIABLES 

c  tauo-sional  decorrelation  time 
c 

COMMON  /PSO/Ct  >  TAUO.FO 

PSOFF»  t .  $44*  TALK)/ <».*$.  §72*  ( FrTAUQ )  **2  »  **2 

return 

END 

SUBROUTINE  RANCOFfXR.Xl.AMP) 

THIS  FUNCTION  RETURNS  a  COMPLEX  RANDOM  NUMBER.  (XR.XU.  WHERE  THE  REAL 
AND  IMAGINARY  PARTS  ARE  INDEPENDENT  ZERO  MEAN  NORMALLY  DISTRIBUTED 
RANDOM  VARIABLES  EACH  WITH  A  VARIANCE  OF  .£*AMF**£. 

arguments 

Xft»RE turned  real  random  sample 

e  It  ^RETURNED  l  MAO 1 NAR  Y  RANDOM  SAMPLE 

C  AMP*SOUAftg  ROOT  OF  TWO  TJMES  THE  VARIANCE  OF  *R  AND  XI 
C  REWIRED  FUNCTIONS 

C  SORT  (FLOATING  POINT  SWARF  ROOT ) 

c  au>>  ^floating  point  natural  loo> 

C  RAND  (SAMPLE  QF  UNIFORM  DISTRIBUTION  BETWEEN  ZERO  AND  0NE> 

C  COS  t FLOAT  I NO  POINT  COSINE  ( 

i  sin  (Floating  point  sine» 

e 

AM-AMP*SCRT  t -ALOG{ RANDf G»  » > 

AN=i . *3 3 IS? *RAND ( 0  * 

XR=AM*cOS t AN) 

Il=AH*SlNtAN> 

RETURN 

END 

FUNCTION  RANDf  t» 

e 

c  Ttt{S  FUNCTION  GENERATES  RANDOM  NUMBERS  UNIFORMLY  DISTRIBUTED 

e  IAPPmOfimatcly)  in  The  interval  of  zero  to  one.  THE  USE*  should 
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1%)  "<*&  (fi>  '<> 


C  REPLACE  THIS  ALGORITHM  WITH  THE  BETTER  ONE (PRESUMABLY)  ON  HIS  OWN 
C  SYSTEM. 

e 

C  ARGUMENTS 
C  I -GUMMY  ARGUMENT 
C  REQUIRED  COMMON  VARIABLES 

C  SEED-RANDOM  NUMBER  S££D< INITIAL  VALUE  MUST  BE  GOO) 

C  REQUIRED  FUNCTIONS 

C  OMGO  (DOUBLE  PRECISION  MOO  FUNCTION) 

C 

DOUBLE  PRECISION  SEED 
COMMON  /RAN/SEED 

SEEP«BM0D(4.  0823 1 25D*07*SEED. 2. 684354360*0$  > 

RAND- <  SEED/2. 684334370*08) 

END 

FUNCTION  PCWER(D) 

C 

C  THIS  FUNCTION  CALCULATES  the  GENERALIZED  POWER  SPECTRUM  AT  A  DELAY  OF 
C  0  INTEGRATED  OVER  ALL  FREQUENCY. 

C 

C  ARGUMENTS 

C  0-DEl.AYCSECl 

C  REQUIRED  COMMON  VARIABLES 

C  Cl -delay  parameter 

C  FO-FRECHJCNCY  SELECTIVE  bandwidth 

C  REWIRED  FUNCTIONS 

e  SORT  (FLOATING  POINT  SQUARE  ROOT) 

C  AB'J  (FLOATINi;-  POINT  ABSOLUTE  VALUE) 

C  EXP  (FLOATING  POINT  NATURAL  EXPONENTIAL) 

e 

COMMON  /PSO/et.TAUO.FO 
FPRlHE*P0*iOftT<  t.*et**2) 

TPFPD-fc. 283-FPR {Mf*D 

|*.70?|*(CI-TPPPD/Ci) 

T»l./(|.*.47047«ABS(Z)) 

PT*t  t .  ?47C«fS6*T-.  0*?3?G3l*T*.  $480242)  *T 
|F(f.LT.O.  )uct  TO  10 

PGWgR=3.  | 4 t t*FFR !»*;*€ XPt  =. S* ( TPFPO/e I  ) **2 ) *FT 
RETURN 

10  P>;»UER-3.  |4l6*FP«tMe*€»PI.5*0l»*2»TPPPt>)*(2.'FT*ExPt-2**i)» 

RETURN 
END 

FUNCTION  TCORtT) 


TlitS  FUNCTION  CALCULATES  ?*€  COMPLEX  SIGNAL  TIME  DECORRELATION 
FUNCTION. 

ARGUMENTS 

T=TJHC  DISPLACEMENT ♦ SEC ) 

RewiPfg  common  variables 

TAUG-S  IGnA?  DE c  OfiRfeL  AT  t  ON  TIME  (EEC  » 

NOrfLAT  FADE  F\AjJ 


S3 


C  NP-I.  FLAT  Ff.JlNO 

C  ND.OT.I,  FREQUENCY  SELECTIVE  FADING 
C  REOU1REB  FUNCTIONS 

C  EXP  (FLOAT I NO  POINT  NATURAt  EXPONENTIAL! 

C  ASS  (FLOAT  INS  POINT  ABSOLUTE  VALUE  I 
C 

COWON  /PSB/Cl.TAUO.FO 
COWON  /OAT/nFP.CFP.NOP.KF>+NO.DS 
C 

C  CHECK  FOR  FLAT  F«OINO 
C 

IF (NO. EQ. I >00  TO  10 
TCOR«£iSP<-<  f/TAU0>**2> 
return 

1C  TC0P»2. 1 46«A5S ( T I  /TAUO 

TCOR*=EXP( -TCORJ  •  U .  •TCOR> 

RETURN 

END 

SUBROUTING  &ETD(TASL.NT» 

c 

C  SUBROUTINES  SETft.  PFT >  AND  POD  ACCOMPLISH  THE  PAST  FOURIER  TRANSFORM 
C  ON  A  COMPLEX  ARRAY. C.  CR  AND  Cl  UILL  b£  USED  AS  A  SHORT  HAND  TO 
C  REPRESENT  THE  REAL  AND  IMAGINARY  PARTS  OF  C.  LET 
C  CR(!>  *  Cd*I-t) 

c  cud  »  cti*n 

C  DEFINE  THE  CORRESPONDENCE  BETWEEN  C.  C«.  AND  Cl.  THE  RESULTS  ARE 
C  RETURNED  IN  THE  c  ARRAY.  THE  NUMBER  OF  POINT?  TO  8>E  TRANSFORMED  IS 
C  NP=;**NT.  THE  FAST  FOURIER  TRANSFORM  RE -.‘U (RES  A  TAfcgfc  REFRiSENTgO  f?Y 

C  THE  COMPLEX  TABL  ARRAY.  MINIMUM  LIMITS  FOR  THE  C  AND  TASL  'FORAYS  ARE 
C 

C  DIMENSION  C(NP*NP»,TA*MNP> 

c 

c  (CR(D».et(D»>  at  the  output  eooals  the  sum  over  i  from  t  to  n?  of 
c 

C  teR(D*C'?S<?.*Pt*M»N^NPV-Cl<D*SI«(2.»Pl««*N^NPi.Cfe(  D*SIN(i.«^*i*M-N/ 
C  HP>*e|(I>*C00(2.*F‘I*M*N.:'NPD 

e 

c  WHERE  THE  COMPLE*  NOTATION  HAS  BEEN  0?FD.  (RgAL.  {MAOIHASV! .  AND 
C  WHERE  N’M-’t  WHEN  J.OE.  t.ANtt.U.l.T.NP/?*! 
e  N^O'-NP-t  WHEN  U.OE  NP/5*I.AND.D.tE.NP 

c  M*t«|  WHEN  I.OE.  I. AND.  t.gT.NP/;'*{ 

£  M*J-NP'*|  WHEN  i.0£.NF'<>2»t.AND.  t.LE.NF- 

e 

C  SUBROUTINE  SEED  INITIALISES  THE  YAfA,  ARRAY  AND  NEED  EE  CALLED  ONLY 
0  WHEN  NT  fS  CitANOCO.  SUBPuuf JN£S  PET  AND  FDD  WHEN  CALLED  IN  lEOUENCS 
f  ACCOMPLISH  THE  TRANSFORM. 

e 

E  CALL  SETDi TA&L.Nf 1 

c 

f  CALX  PFT(C.TAEU‘ 

e  call  RDDtei 

c 


St 


C  IN  THE  ABOVE  SUM.  M  IS  PROPORTIONAL  TO  THE  VARIABLE  REPRESENTED  BY  THE 
C  I  INDEX.  FOR  EXAMPLE.  IF  (CR< I ) .01 < I ) >  IS  A  FREQUENCY  COEFFICIENT. 

C  THEN  ZERO  FREQUENCY  IS  REPRESENTED  BY  M  EQUALTO  ZERO  AND  I  EQUAL  TO 
C  ONE.  SIMILARLY.  N  IS  PROPORTIONAL  TO  THE  VARIABLE  REPRESENTED  BY  THE 
C  U  INDEX. 

C 

C  FAST  FOURIER  TRANSFORMS  ASSUME  CONTINUITY  OF  THE  FUNCTION  TO  BE 
C  TRANSFORMED  AND  THE  RESULT  AT  I<OR  J>  EQUAL  TO  NP/2+1  WHICH  REPRESENTS 

C  THE  EXTREMES  OF  M(OR  N).  THE  FUNCTIONS  ARE  ALSO  CONTINUOUS  BETWEEN 

C  <CR( 1 >  >CI ( 1 ) )  AND  ( CR  <  NP ) .  C I  ( NP ) ) . 

C 

C  REQUIRED  FUNCTIONS 

C  FLOAT  (FIXED  TO  FLOATING  POINT  CONVERSION) 

C  MOD  (FIXED  POINT  MOD) 

C  COS  (FLOATING  POINT  COSINE) 

C  SIN  (FLOATING  POINT  SINE) 

C 

C  #*#**###*«******•*#«***•*»*»*«*'.'**»««****»•**••*****************•***** 

C  THE  ARRAYS  IN  SETD.PFT,  AND  PDD  HAVE  BEEN  DIMENSIONED  TO  ONE 
C  INTERNALLY.  THIS  MAY  NOT  BE  ALLOWED  ON  THE  USER  SYSTEM.  IT  MAY  BE 
C  NECESSARY  TO  DIMENSION  TO  THE  LIMITS  IN  THE  MAIN  PROGRAM  OR  TO  USE 
C  VARIABLE  ARRAY  DECLARATORS. 

C  •*••**»•»»***«««****•»••••*«••*»•*****•»***»*•****•»•****•**••*•*«»«»» 

c 

DIMENSION  TABL<1) 

COMMON  /FFT/N2-NTI.NL 
DATA  PI/3,1415926536/ 

NTI=NT 

N2»2**NT 

NL=N2-1 

NN=N2/2-l 

TP I =2.  #P  I  /FLOAT  (N2) 

TABU  1  )  =  1 . 

TABL(2)»0. 

DO  40  1=1. NN 
J»I  +  I 
11*0 

DO  20  N=1 » NT 
11*11*11 >MOD  <  J 1 2 ) 

20  J*U/2 

K=I  +  IU 

TABL  <  K ) =COS ( TP I *FLOAT  <II)> 

40  TABL(k>1 )*SJN(TPI*FLOAT( II) ) 

RETURN 

END 

SUBROUTINE  PFT(C.TABL) 

DIMENSION  C(l).TABLd) 

COMMON  /FFT/N2P.NTI.NL 

Nl*2 

N2»N2P 

DO  1  I«1.NTI 

NN»-N2-N2 


-T))S 


DO  2  J=2,N1,2 
NN=NN+N2+N2 
N3=NN+N2 
DO  3  K=2> N2* 2 
N4=N3+K 

STR=C ( N4- 1 >  *TABL ( J- 1) -C ( N4 ) *TABL ( J ) 
STI=C(N4-1 )*TABUJ)+C<N4)#TABL<U-i) 
N5=NN+K 

C(N4-1 )=C(N5-1 >-STR 
C(N4)=C(N5) -ST I 
C(N5-1 >=C(N5-1 )+STR 
3  C(N5)=C<N5)+STI 

2  CONTINUE 
N1=N1+N1 

1  N2=N2/2 

RETURN 
END 

SUBROUTINE  PDD(C) 

C 

C  REQUIRED  FUNCTIONS 
C  MOD  (FIXED  POINT  MOD) 

C 

DIMENSION  C < 1 > 

COMMON  /FFT/N2.NTI.NL 
DO  l  1=1, NL 
J=I 
11=0 

DO  3  N=1,NTI 
11=11+1 I*MOD( J« 2 ) 

3  J-J/2 

IF( 1 1  —  I  )  1, 1 ,7 
7  IN-I+I+1 

IIN-II+II+l 
OR=C( IN) 

DI=C( IN+1 > 

C( IN )=C  < I IN ) 

C< IN+I )=C< I IN+1 ) 

C(IIN)=DR 
CUIN+I  )=DI 
1  CONTINUE 

RETURN 
END 
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